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We measured the absolute lengths of three single crystal silicon samples by means of an imaging 
Twyman-Green interferometer in the temperature range from 7 K to 293 K with uncertainties 
of about 1 nm. From these measurements we extract the coefficient of thermal expansion with 
uncertainties in the order of 1 x 10~^/K. To access the functional dependence of the length on 
the temperature usually polynomials are fitted to the data. Instead we use a physically motivated 
model equation with 7 fit parameters for the whole temperature range. The coefficient of thermal 
expansion is obtained from the derivative of the best fit. The measurements conducted in 2012 and 
2014 demonstrate a high reproducibility and the agreement of two independently produced samples 
supports single crystal silicon as reference material for thermal expansion. Although the results for 
all three samples agree with each other and with measurements performed at other institutes, they 
significantly differ from the recommended values for thermal expansion of crystalline silicon. 


I. INTRODUCTION 


High accuracy knowledge of the coefficient of thermal 
expansion (GTE) of material samples is essential for de¬ 
velopment and characterization of ultra stable materials 
needed e.g. in semiconductor industry, precision optics 
or in aerospace applications^’^. To enable high accuracy 
operation many measurement systems use reference ma¬ 
terial samples as calibration standards. For this purpose 
single crystal silicon (SGS) is a commonly used material 
that is particularly suited. SGS is available as off-the- 
shelf product with high purity. Together with its sin¬ 
gle crystal structure this guaranties supplier independent 
thermal expansion, that is isotropic and low compared to 
metals. Further it offers a high thermal conductivity, es¬ 
pecially towards cryogenic temperatures, which supports 
a homogeneous temperature distribution over the sam¬ 
ple. 

Thermal expansion of ultra stable materials at cryo¬ 
genic temperatures is of increasing scientific interest 
driven by technological applications^^^. High accuracy 
GTE values at cryogenic temperatures are essential e.g. 
for construction and operation of spaceborn telescopes 
as the Herschel Space Observatory® (operated at 85 K), 
the James Webb Space Telescope® (to be operated be¬ 
low 50 K), or the SPIGA telescope^® (desired: 5 K). In a 
round robin performed a few years ago under the guid¬ 
ance of the European Space Agency (ESA) and the Gen- 
tre Spatial de Liege (GSL) the thermo-mechanical char¬ 
acterization abilities of some of the world-wide best in¬ 
stitutes were compared^ A clear lack of high accuracy 
measurement capabilities at cryogenic temperatures was 
pointed out. For the design of ultra stable structures 
it is necessary to characterize the thermal expansion of 
these materials with an uncertainty below 3 x 10“®/K 
in the entire in-service temperature range^®. Thus, for 
a reference material as SGS a similar or even lower un¬ 
certainty is indispensable. To realize this PTB’s Ultra 
Precision Interferometer (UPI) was equipped with an ex¬ 
tended measurement pathway, enabling absolute length 
measurements from 7 K to 293 K®. 


A review of measurement techniques for the GTE of 
metals and alloys is given by James et al.®®. Often the 
GTE is determined from the ratio of length change over 
temperature change Al/AT (as e.g. by Lyon et al.®^) re¬ 
quiring additional corrections to account for the finite 
non-zero interval size®®. Instead of measuring length 
changes, we measure the absolute length as a function 
of temperature. Further, rather than fitting polynomials 
in several regions to the data, we use a physically moti¬ 
vated model equation, that covers the whole investigated 
temperature range, with only seven free parameters. By 
this we avoid an overestimation of the data, and give a 
coherent description of the functional dependence. 

In this paper we report on GTE measurements on three 
different samples, manufactured from two independent 
suppliers. All measurements agree with each other, and 
with other high accuracy measurements. Albeit a com¬ 
parison with the GTE values recommended for crystalline 
silicon unveils significant systematic deviations of 3 to 4 cr 
in a wide temperature range (90 K to 210 K). 


II. EXPERIMENTAL SETUP 


The experimental setup has been described in detail 
earlier®. For completeness we recall the main aspects 
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FIG. 1. Schematic of the beam path in PTB’s Ultra Precision 
Interferometer. Prismatic samples are placed in the extended 
probe beam on the left. Its temperature can be varied between 
7 K and 293 K by a pulse tube cryocooler. 
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in short. The measurements were performed with PTB’s 
UPI, an imaging Twyman-Green interferometer operated 
under vacuum (Fig. 1). The expanded beam from three 
alternatively used frequency stabilized lasers is split up 
into a reference and a probe beam. The probe beam is 
directed towards the sample, and the arising interference 
pattern is imaged on a CCD-Array. The length measure¬ 
ments are performed by phase stepping interferometry. 
For each length measurement 10 pictures of the fringe¬ 
less interference pattern at 10 different phase steps are 
recorded. The length is determined twice, from 5 pic¬ 
tures for each value, and both values are averaged. Phase 
stepping is realized by tilting the compensation plate in 
the reference beam by a piezo stepper. 

For measurements at cryogenic temperatures the probe 
beam is elongated by an extended measurement pathway 
(EMP), housed in a separated vacuum chamber. The 
cooling is achieved by a pulse tube cryocooler (PTC) 
with two cooling stages, enabling the variation of the 
sample temperature from 7 K to 293 K. The samples 
are surrounded by a copper shield and measurements are 
performed under a low pressure (p « 1 mbar) helium 
atmosphere that guaranties a homogeneous temperature 
distribution. The optical path dilatation due to the re¬ 
fractive index of diluted helium is regarded as described 
earlier^. The temperature is measured by means of two 
Rhodium-Iron-Resistance-Thermometers (RIRT), which 
are fitted into holes in the center of the samples and 
dummy-samples respectively. The RIRT sensors have 
been calibrated by PTB and their measurement uncer¬ 
tainty rises from 15 mK at a few K to 25 mK at 293 K. 

At temperatures below 35 K the PTC is still oper¬ 
ating during measurements. To minimize the influence 
of vibrations the camera is triggered by the frequency 
of the PTC. At temperatures above 35 K the PTC is 
switched off and measurements are performed while the 
temperature slowly drifts upwards. This restricts the size 
of temperature intervals that can be investigated within 
one day. But since the absolute length of the samples 
is measured, it is not necessary to perform all measure¬ 
ments in one run. Indeed within one cooling cycle an 
arbitrary chosen temperature interval of about 15 K was 
investigated. By this also the emergence of deposits on 
the surfaces of samples and base plate, that were ob¬ 
served at longer times (> 1 day) at low temperatures, 
was prevented. 


III. MEASUREMENTS 

We have performed two measurement series: one in 
2012 on sample #1 and one in 2014 on samples #2 
and #3. Samples #1 and #2 have been prepared from 
SCS material which was acquired from Wacker Siltronic, 
Germany by PTB some years ago, as mentioned by 
Becker^® and Schodeh^. These samples were cut in (100) 
orientation from the 0-Zone of a high purity dislocation 
free float-zone(FZ)-silicon crystal. The level of impurities 


was found very low (oxygen: (1... 2) x 10^® cm“^, car¬ 
bon: (2... 6) X 10^® cm“^, nitrogen: < lO^"* cm“^) and 
the 0-Zone is free from extended swirl defects. Further 
information on the material can be found in^®. The gauge 
block shaped samples have a length of about 35 mm and 
two parallel faces of 9 mm x 20 mm cross section. The 
latter were lapped by Kolb & Baumann GmbH & Co. 
KG (KoBa) to optical quality. Sample #3 was acquired 
by PTB together with a “Reference Material Certificate”, 
containing certified CTE values for this sample, from the 
National Metrology Institute of Japan (NMIJ)^®. This 
sample was as well made from a high-purity FZ-silicon 
crystal produced by Shin-Etsu Handotai Co., Ltd. This 
third sample is a rectangular block with a length of 
30 mm and two parallel faces of 10 mm x 10 mm (Fig. 
2). These faces have been lapped in the same way as for 
the other two samples at KoBa to optical quality. For 
the measurements the samples were wrung onto a base 
plate also made from SCS and polished in the same way 
by KoBa. 

In the 2012 measurement series not always optimal 
conditions have been used as was discovered later on. 
Reasons for suboptimal conditions were: Deposits from 
condensation due to the presence of residual air or small 
leakages in the vacuum system caused by temperature 
gradients. Measurements that were subject to deposits 
were suspended from the 2012 measurement series. Also 
a few measurements had to be suspended due to high lo¬ 
cal or temporal temperature gradients, that were caused 
by unilateral or too strong heating. 

In the 2014 measurement series high local or tempo¬ 
ral temperature gradients were avoided and deposits from 
condensation were prevented by using a protective helium 
atmosphere during cool down and between measurements 
as well as a cold trap operated at 77 K with liquid nitro¬ 
gen. The protective helium atmosphere first of all reduces 
the pressure difference to ambient air pressure and by 
this also leakage amounts. Second of all it dilutes resid¬ 
ual gases and prevents their condensation. Additionally 
the pressure gradient between the EMP inner vacuum 
chamber and the surrounding isolating vacuum chamber 
is inverted, and leakages from the surrounding vacuum 
chamber are suppressed completely. For the pressure 
of the protective helium atmosphere a high pressure of 
900 mbar turned out to be most effective. A side effect 
was the shorter time needed to cool down, due to better 
thermal linking by the heat transport of the protective 
atmosphere. In general the times at cryogenic tempera¬ 
tures were kept short, i.e. typically under 8 hours and the 
protective helium gas was exchanged about every three 
hours, to dilute contaminations. 

The temperature sensors described in Sec. II were 
placed in drillings in the center of sample #1 and #2. 
Sample ^3 is too small to contain a temperature sensor 
thus its temperature is determined as average of two sen¬ 
sors placed in two silicon blocks that were wrung to the 
same base plate (one of them is sample #2). Measure¬ 
ment uncertainties are discussed in detail in Appendix A. 
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IV. CTE DETERMINATION 


The coefficient of thermal expansion is determined ac¬ 
cording to the ISO definition^® 


a{T) = 


1 dl{T) 
Ikt dT 


( 1 ) 


with the temperature dependent length l(T) and the 
length Zrt at room temperature RT = 293.15 K. 

To determine the CTE according to Equation 1, the 
functional dependence of the length on the temperature 
needs to be known. This can be achieved by fitting an 
appropriate function to the data. Thereto usually poly¬ 
nomials are used and fitted in several regions to the data. 
But of course selecting polynomial order and fitting re¬ 
gions also introduces uncertainties. Thus to obtain a pos¬ 
sibly low uncertainty, a functional description over the 
whole temperature range, that is derived from physics is 
desirable. 

In a simple model, that neglects phonon dispersion, 
the CTE can be expressed according to the Griineisen 
equation®® 


a(^) = ^KjCy(T), 


( 2 ) 


as being proportional to the specific heat c„(T) per unit 
volume at constant volume. Here also the temperature 
dependencies of the Griineisen parameter 7 and harmonic 
compressibility k are neglected. A simple model for the 
temperature dependence of the specific heat is provided 
by the Einstein model of specific heat: 

rrr.^ 

Cv{ ) y j (_i^geE/T)2- (3) 


Here Qe is the material dependent Einstein- 
Temperature. Equations 2 and 3 represent only an 
approximate functional dependence of the GTE on 
the temperature. But they can serve to constitute an 
appropriate fit function by summation over several so 
called “Einstein-Terms”®®, such that the CTE can be 
described by 


.. m 


ak 




\T J (-l-ke»'=/'^)®’ 


( 4 ) 


and the fit function is obtained by integration to be 


1{T) — Iq + ak 


fc=i 


— I -|- 


( 5 ) 


The fit parameters 0fc have the unit of temperature but 
do not directly represent physical properties, although 
they are roughly related to peaks in the frequency spec¬ 
trum of the phonon density of states. The fit parameter 
lo corresponds to the sample’s length at T = 0 K. For 
the measurements considered here on SCS to = 3, i.e. 
n = 7 fit parameters are sufficient to fit the data very 
well within the expected uncertainty. Fitting and un¬ 
certainty evaluation of the CTE values are described in 
Appendix B. 
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FIG. 2. a) Measured length versus temperature for sam¬ 
ple #2, fit (red line) according to Eq. 5, and below residuals 
from the ht and overall uncertainty of the length versus tem¬ 
perature data (Appendix A), b) Deduced CTE according to 
Eq. 4 and overall uncertainty (Appendix B). 


V. RESULTS 


The length versus temperature data obtained for sam¬ 
ple #2 are plotted in Fig 2a together with the fit (Equa¬ 
tion 5) and its residuals. The latter are small compared 
to the overall uncertainty (red) estimated in Appendix A, 
demonstrating the appropriateness of the fit function. 
The fit parameters are tabulated in Tab. I in Appendix B 
and the results obtained for sample #2 are tabulated in 
Tab. H in Appendix C 

provided online as supplemental material. The corre¬ 
sponding CTE values are plotted in Fig 2b (Equation 4) 
and below its overall uncertainty (blue), whose estima¬ 
tion is described in Appendix B. 

The results for all three samplesare compared in Fig 3. 
Plotted is the respective CTE as a band (half width of 
1 a), minus the average of all three measurements. All 
three measurements agree very well within their uncer- 
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FIG. 3. Comparison of all three measurements performed 
at PTB: In 2012 on sample (light blue), and in 2014 on 
samples (dark blue) and ^3 (purple). Plotted is a band 
with a half width of 1 ct around the respective GTE minus 
the average GTE of all three measurements. 

tainties. The largest difference of two of the measure¬ 
ments is smaller than 3 x 10“®/K. The best results are 
obtained for sample #2, since it was investigated with the 
optimized measurement procedure. Furthermore in com¬ 
parison to sample ^3, that was investigated in the same 
measurement series, it contained a temperature sensor 
placed in its center, whereas the temperature of sam¬ 
ple #3 was determined as average of two sensors placed 
in adjacent samples, as described in Sec. III. 

A comparison with other^^ high accuracy GTE mea¬ 
surements performed on crystalline silicon is presented 
in Fig 4. Plotted is a band with a half width of 1 tr 
around the respective GTE, minus the GTE values rec¬ 
ommended by GODATA as reference data^^. The results 
from this work (sample #2) are plotted in blue, the data 
from the reference certificate acquired from NMIJ^® in 
red and measurements from the Jet Propulsion Labora¬ 
tory (JPL) in green®. The l-cr-uncertainty of the refer¬ 
ence data is indicated by the gray shaded band accord¬ 
ing to the original publication^^, that is commonly cited. 
Although it is worth to be mentioned, that the data in 
the recommendation of GODATA^®, are provided with 
an expanded uncertainty of 2 x 10“®/K for temperatures 
between 40 K to 300 K, which is up to 20 times larger 
than the uncertainty stated originally. 

All but the recommended reference data agree within 
their uncertainties in the whole temperature range. But 
in the temperature range from 90 K to 230 K the recom¬ 
mended reference data differ from the results of this work 
by about 4 x 10“®/K or 3cr to 4a. Here it is noticed, that 
the measurements^'* which are the basis for the reference 
data at temperatures below 300 K, were performed on 
a polycrystalline but not monocrystalline silicon sample, 
as in all other measurements presented here. 

The temperature measurements in this work refer to 
the ITS-90^^, while for the recommended reference data 


FIG. 4. Gomparison of high accuracy GTE measurements 
performed on crystalline silicon. Plotted are bands with a 
half width of 1 cr around the respective GTE minus recom¬ 
mended reference data^^. Gompared are: The results from 
this work on sample ^2 (blue), the data from the reference 
certificate acquired from NMIJ*® (red), measurements from 
JPL® (green), and the l-u-uncertainty of the reference data 
(gray). 

the IPTS-68 was used, as mentioned by Kroeger and 
Swenson*®. Taking this into account reduces the mis¬ 
match of both measurements slightly, but not signifi¬ 
cantly. The differences between thermodynamic temper¬ 
atures and the ITS-90 were estimated by Fischer et al.^® 
and are smaller than 9 mK in the temperature range 
considered in this work, which is less than half the tem¬ 
perature uncertainty (Sec. II) and thus negligible for our 
results. 

The significance of the mismatch between the recom¬ 
mended reference data and this work becomes more ev¬ 
ident considering the thermal strain s(T) = (l(T) — 
Irt)/Irt- In Figure 5 the recommended thermal strain 
data are compared with the results from this work. While 
the strain difference of all three measurements from this 
work is compatible with the estimated uncertainties, the 
strain difference between this work and the recommended 
data increases towards lower temperatures and becomes 
as large as 13 a, referring to the measurement on sam¬ 
ple #2. As the strain is proportional to the integration 
of the GTE, the difference visible in Fig. 4 is now accu¬ 
mulated, emphasizing the significance of the discrepancy. 
This also illustrates the advantage of absolute length 
measurements, since the lengths measured at different 
temperatures can be directly compared and do not re¬ 
quire an integration. 

VI. CONCLUSION 

The GTE of three SGS samples were determined from 
absolute length measurements at temperatures between 
7 K and 293 K. Uncertainties in the order of or smaller 
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FIG. 5. Comparison of thermal strain results. Plotted are 
bands with a half width of 1 cr around the respective strain 
minus the strain measured for sample #2. Compared are: 
The results from this work on sample (light blue), on 
sample ^2 (dark blue) and on sample #3 (purple) with the 
recommended reference data^'* (black line), which are pro¬ 
vided without uncertainty. 


than 3 x 10“®/K were obtained (10 K to 290 K). The 
two independent measurement series performed in 2012 
and in 2014 demonstrate a high degree of coincidence, 
no matter if CTE or thermal stain are compared. Fur¬ 
thermore one of the two samples investigated in 2014 was 
manufactured in Japan, completely independent from the 
other two samples manufactured in Germany. Thus also 
the universality of SCS as reference material is demon¬ 
strated with an accuracy of a few times 10“®/K. 

A comparison with other high accuracy measurements 
performed in the last decade demonstrates a good agree¬ 
ment. Albeit a mismatch is found with the reference 
data recommended by CODATA^^’^^ in a wide tempera¬ 
ture range (90 K to 230 K). The CTE difference in this 
region corresponds to 4 x 10“®/K, which is about four 
times the uncertainty. This mismatch is supported by 
the measurements at the other institutes. In recent years 
accuracies below 3 x 10“®/K are coming more and more 
to the fore in industrial applications, thus reliable refer¬ 
ence data are mandatory. As SCS is commonly used as 
reference material for thermal expansion, also the recom¬ 
mended reference data should be based on SCS. 
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Appendix A: Measurement uncertainty 

The uncertainty of the individual data points stems 
from both, the length and the temperature measure¬ 
ments. In general each of these contributions consists of 
a type A (evaluated by statistical methods) and a type 
B (not evaluated by statistical methods) part^®. The un¬ 
certainty of the pure length measurements is evaluated 
as described in^. The type A uncertainty is estimated 
to be ma( 0 = 1 nm and the type B uncertainty to be 
wb( 0 = 0-42 nm. In the case of sample #2 a settling of 
the wringing contact of 0.039 nm/day was observed, as 
apparent from length measurements performed at room 
temperature. The settling was corrected for. 

The uncertainty of the temperature measurements is 
treated as type B uncertainty. It is estimated according 
to the GUM-rules^® for a rectangular distribution from 
the temperature difference of both sensors and the sensor 
uncertainty to be 

mb(T) = + (0-016+ 3.3-10-5 r)A (Al) 

The uncertainty of temperature u{T) is propagated to 
the length dimension by 

= “B(r)\/(J:)“ + «(|:) . (A2) 

with u(^) = 2 X 10“®/rt being the estimated uncer¬ 
tainty of the derivative to account for the uncertainty 
that remains when the derivative equals zero. 

The total uncertainty of the length versus temperature 
is obtained from the sum of squares to be 

u{l) = yj Ua(0^ + Mb(0^ + UB,Tiiy, (A3) 
and is plotted in Figure 2a for sample #2. 


Appendix B: CTE uncertainty 

Equation 5 is fitted to the length versus tempera¬ 
ture data using the Levenberg-Marquardt algorithm with 
Mathematica 8.5 (Wolfram Research) by the built-in 
function NonlinearModelFit. For fitting, the data are 
weighted according to the combined measurement uncer¬ 
tainty from the individual length and temperature mea¬ 
surements with l/u^ = 1 /('ua(0^ + wb,;(0^ + wb,t(0 ^- 

Some of the n = 7 fit parameters pi (Tab. I) are 
strongly correlated. Hence, to estimate the type A un¬ 
certainty of the fit the law of propagation of uncer¬ 
tainty, for correlated input quantities has to be applied^®. 
This requires, to take into account both, the variances 
Cii = u{pi)'^, and also the covariances Cy of the param¬ 
eters Pi and pj, being the diagonal and the off-diagonal 







6 


TABLE I. Values of the fit parameters obtained for sample #2 
and corresponding statistical uncertainties. 


Parameter 

Value 

u{pi) 

Pi 

= Ml 

-3.398 X 10"'^* 

6.3 X 10”"" 

P2 


199.61 

0.20 

P3 

= 02 

1.487 X lO”®^ 

3.2 X 10”^° 

P4 

= 62 

612.00 

0.61 

P5 

= 03 

3.496 X 10“°* 

4.4 X 10““ 

P6 

= 6^ 

890.05 

0.09 

P7 

= h 

0.0348286997 

1.5 X 10““ 


elements of the n x n covariance matrix of the fit. Thus, 
the type A uncertainty of the fit is given by 


ua{KT)) 


n n 


EE 


dl dl 


(Bl) 


in which dl/dpi are the partial derivatives of the Fit 
(Eq. 5) with respect to the parameter pi. Analogously 
the type A uncertainty of the CTE is given by 


ua(q;(T)) 




EE 


da da 


(B2) 


with dajdpi being the partial derivatives of Eq. 4 with 
respect to the parameter pi. 

Additionally a type B error can occur due to inappro¬ 
priateness of the fit function. To analyze this we consider 
the adjacent median of the fit residuals. From its com¬ 
mon mode amplitude we deduce an estimate of the type 
B error of mb .modei(^T)) = 1 nm. Since the GUM^® does 
not offer a simple model for the propagation of the type 
B uncertainty of a quantity to its derivative, the type B 
uncertainty of the fitted length 

Ub(K^)) = Mb(0^ + Wb.t( 0^ + WB.model(^(T))^ (B3) 


At last we consider a circumstance inherently con¬ 
nected with the determination of the slope of measured 
data. The closer two data points are, the harder it is 
to tell the slope of a line connecting them. As the line 
could go from the lower end of the uncertainty inter¬ 
val of one point, to the upper end of the uncertainty 
interval of the other, and vice versa. Since the mea¬ 
sured data are distributed over a wide temperature range, 
this becomes only crucial at the borders of the investi¬ 
gated temperature range. Especially for low tempera¬ 
tures, when T approaches zero, according to equation 4 
applies lim a{T) = 0. As a matter of principle the 

T —^0 

same applies for the type A uncertainty (equation B2): 
^ini UA{a) = 0 and the type B uncertainty (equation B4), 

i.e. ^im mb (a) = 0. This is a property of the fit function 

used, which is only an approximation and reality might 
differ. Although it can be derived from the third law of 
thermodynamics, that a{T = 0) = 0, for small tempera¬ 
tures above T = 0 the true shape might differ from the 
fit function. Here it is noticed, that in this temperature 
region small positive CTE values as large as 1 x 10“®/K 
at T = 14 K have been reported^^. As this is in the 
order of the uncertainty, it cannot be validated by the 
measurements presented in this work. We estimate this 
slope uncertainty to be 


Ud{a{T)) 


V2ua{1) 

Zrt(AT-2 T-T)’ 


(B5) 


here ma( 0 is the type A uncertainty of the individ¬ 
ual length measurements, AT = Tmax — Tmin is the 
size of the temperature range under investigation and 
T = Tmin + at/ 2 is the central temperature. The slope 
uncertainty Ud{a{T)) is small in the largest part of the 
investigated temperature range, but becomes the dom¬ 
inating contribution below 10 K and above 290 K and 
avoiding an overestimation of the data in these regions. 

All three contributions are combined to give the overall 
uncertainty of the CTE 


u{a) = a/ma(q;(T))2 -|- mb(q;(T))2 -|- ud{a{T)Y , (B6) 


is projected via the ratio of the type A uncertainties of 
the CTE and the length (equations B2 and Bl), such 
that the correlation is regarded here as well 

mb(«(T))= ^^^^^^E b(Z(T)). (B4) 

ma(/(T)) 


that is plotted in Figure 2b for sample #2. 


Appendix C: Tabulated CTE values 
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( 2011 ). 

® S. Stoupin and Y. Shvyd’ko, Phys. Rev. Lett. 104 , 085901 



7 


TABLE II. Tabulated CTE values a /(10“®/K) and uncer¬ 
tainties u(a) /(10~®/K) versus temperature T /K, extracted 
from the measurements on sample #2. 


T 

a 

u(a) 

8.15 

0.0 

11.2 

13.15 

-0.1 

3.0 

18.15 

-2.0 

1.7 

23.15 

-13.1 

1.3 

28.15 

-40.9 

1.4 

33.15 

-86.2 

2.0 

38.15 

-144.0 

2.5 

43.15 

-207.9 

2.4 

48.15 

-271.9 

2.0 

53.15 

-331.7 

1.5 

58.15 

-383.7 

1.2 

63.15 

-425.5 

1.3 

68.15 

-455.5 

1.5 

73.15 

-472.6 

1.6 

78.15 

-476.2 

1.6 

83.15 

-466.5 

1.5 

88.15 

-443.8 

1.3 

93.15 

-409.0 

1.2 

98.15 

-363.0 

1.1 

103.15 

-306.9 

1.1 

108.15 

-242.1 

1.2 

113.15 

-169.6 

1.3 

118.15 

-90.8 

1.3 

123.15 

-6.8 

1.4 

128.15 

81.4 

1.3 

133.15 

172.8 

1.2 

138.15 

266.4 

1.1 

143.15 

361.7 

1.0 

148.15 

457.7 

0.9 

153.15 

554.0 

0.8 

158.15 

650.1 

0.8 

163.15 

745.5 

0.9 

168.15 

839.9 

0.9 

173.15 

932.9 

1.1 

178.15 

1024.4 

1.2 

183.15 

1114.1 

1.3 

188.15 

1201.9 

1.3 

193.15 

1287.6 

1.4 

198.15 

1371.2 

1.4 

203.15 
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